Measuring entanglement growth in quench dynamics of bosons in an optical lattice 
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We discuss a scheme to measure the many-body entanglement growth during quench dynamics 
with bosonic atoms in optical lattices. By making use of a ID or 2D setup in which two copies of 
the same state are prepared, we show how arbitrary order Renyi entropies can be extracted using 
tunnel-coupling between the copies and measurement of the parity of on-site occupation numbers, as 
has been performed in recent experiments. We illustrate these ideas for a Superfluid-Mott insulator 
quench in the Bose-Hubbard model, and also for hard-core bosons, and show that the scheme is 
robust against imperfections in the measurements. 



PACS numbers: 37.10.Jk, 03.65.Ud, 03.67.Mn, 

Entanglement is a basic feature of many body quantum 
systems pQ, and underlies the complexity of simulating 
quantum physics on a classical computer [2]. The expo- 
nential scaling of resources to represent and propagate a 
general many body quantum state on a classical device 
has motivated the development of quantum simulators 
[3], and significant progress has been made in building 
both analog and digital quantum simulators with cold 
atoms and ions for equilibrium and non-equilibrium dy- 
namics. This is exemplified by quantitative measurement 
of phase diagrams, studies of quantum phase transitions, 
and quench dynamics. An outstanding challenge, how- 
ever, is direct measurement of (potentially large scale) en- 
tanglement, and monitoring entanglement growth in non- 
equilibrium dynamics. Below we address these questions 
by discussing measurement scenarios for entanglement 
entropies, using multiple copies of a quantum system and 
measurements with a quantum gas microscope [4, 5 . We 
illustrate these ideas in the context of quench dynamics of 
bosons in ID optical lattices. This example is motivated 
by recent experiments [6 , where quench dynamics were 
observable for times not accessible to (classical) t-DMRG 
simulations of Hubbard dynamics [7HIU] due to entan- 
glement growth [TTHT4] . Here the measurement protocol 
will directly reveal this entanglement growth, and simul- 
taneously monitor the purity of the total system state. By 
comparing copies, these tools will also provide a protocol 
for the verification of a quantum simulator. 

We are interested in quantum dynamics of an (ideally) 
isolated quantum system as represented by our atomic 
quantum simulator. In particular, we study a system 
where we prepare an initial state |\I/(0)), which evolves 
with a Hamiltonian # as \V(t)) = exp(-iHt / h)\V (0)) . If 
the system can be divided into two subsystems A and B 
and is in a pure state at time t, p — then the en- 

tanglement of the system can be characterized in terms of 
the entropy of the reduced density matrix, pA = tr#{p}. 
This is commonly computed as the von Neumann en- 
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FIG. 1: (a) Measurement of n — 2 Renyi entropy for bosons 
in an optical lattice. First, two instances of the many-body 
state are produced (shown here for a single site in a ID chain 
or 2D plane). Then tunnelling is switched off within each 
copy, and the barrier between the copies lowered to realize a 
beamsplitter operation between the copies. Finally, the parity 
of the atom number is measured at each site. This measure- 
ment is repeated to obtain expectation values for the swap 
operator V2, from which the Renyi entropy can be computed 
(see text), (b) Example measurement outcome for a single 
shot on a quantum chain. Here the measurement result for 
the whole system swap operator V<£ is 1, since the total 

number of particles in Copy 2 is even. For the swap of the 
first three sites V^ 1 ' 2 ' 3 ^ is —1, since this number is odd. 



tropy Svn(p) = — tr{plogp}. If A and B are in a prod- 
uct state = \^a) ® \^b): then pA will also represent 
a pure state with ti{p 2 A } = 1, and Svn(pa) = 0. Be- 
low we will discuss measurement of the Renyi entropy 
of order n > 2, S n (p) = l°g tr{p n }, which gives 
quantitative bounds for a variety of measures of entan- 
glement, e.g., the concurrence [T5HT8] , and also Svn(p)- 
For example, we know that Svn(p) [= lim Sn(p)], and 

n— vl 

as dS n (p)/dn < 0, S V n(p) > S 2 (p)- From d 2 S n /dn 2 > 
a stronger bound Svn(p) > 262 (p) — Ss(p) can be ob- 
tained [19] if we measure multiple Renyi entropies. 
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In order to measure the Renyi entropy S n (p) for a 
state p, we require n copies of the state prepared in 
parallel, and the possibility to implement operations 
that exchange the n copies [T5l [201424] . As shown in 
Ref. [20 , the quantity tr{p n } can be written [20j [25] 
as tr{p n } = ti{V n p® n }, where the shift operator on n 
copies, Vn^x) . . . \i/> n ) = \^ n )\^i) • • • l^n-i)- Therefore, 
the measurement of the Renyi entropy can be reduced to 
determining the expectation value (V n ) on the n copies. 
Measurements of inner products can be made in this way 
by entangling a state with auxilliary qubits or a quan- 
tum switch [20[ [26H28] , as has been demonstrated for a 
few entangled photons [29H3T] . In Ref. [22 it was shown 
that a beamsplitter operation on two copies is sufficient 
to measure the purity tr{p 2 } for bosonic systems. Our 
study of entanglement growth is based on generalizing 
these techniques to measure Renyi entropies of arbitrary 
order n. Remarkably, the corresponding experimental 
tools (controlled tunnelling between multiple copies and 
measurement of onsite occupation numbers modulo n) 
are now available with single-site addressing in a quan- 
tum gas microscope [4j [5] . 

Below we discuss the protocol to measure S n (p) for 
arbitrary n, and summarize the details of the procedure 
for the simplest example n = 2. We also provide a sim- 
plified scheme for hard-core bosons, and discuss the ro- 
bustness of measurements with respect to imperfections. 
We note throughout that this scheme allows simultane- 
ous measurement of S n (p) for the whole state and for 
every reduced subsystem (i.e., p = Pa) [48 . Thus, while 
measurement of the entanglement in this form is based on 
the assumptions (i) that the initial states for the whole 
system are pure, and (ii) that each evolves under the 
same Hamiltonian, these assumptions can be checked 
directly by measuring S n (p) for the whole system. In 
quench dynamics, these assumptions should be well ful- 
filled in experiments beginning from low-entropy states 
[32] [33] . Moreover, by monitoring the copies over time 
and measuring, e.g., tr{pip2} [20j[25], this scheme pro- 
vides a means to verify a quantum simulator, determin- 
ing whether the evolution of the copies is coherent and 
identical on the level of many-body wavefunctions. 

Scheme for arbitrary n — The procedure to measure 
the Renyi entropy S n (p) consists of three steps: (1) n 
identical instances of the many-body state are prepared 
in parallel (either in n ID chains, or in n planes in 2D). 
This can be performed, e.g., by beginning from a low 
entropy initial state such as a Mott insulator [32j [33] . 
and manipulating the lattice potential identically for the 
two copies, i.e., allowing them to evolve under the same 
Hamiltonian. In this step, the lattice depth between the 
copies must remain large so that these are isolated from 
each other. (2) We then make the lattice deep within 
each copy of the state, to prevent tunnelling, and perform 
a discrete Fourier transform U^ T operation on the copies. 
If the bosonic annihilation operator for site i in copy 



c G {1, . . . , n} is a^ c , then 

i n 

Ur-.aj^^a^-^-V. (1) 

V 1 = 1 

This can be achieved by a successive application of tun- 
nelling between the copies, and shifting the relative po- 
tential depths (to produce elements analogous to beam- 
splitters and phase-shifters in the optical implementation 
of this operation [34 J. This operation is very simple for 
small n, as discussed below, and should be performed 
on the whole system in parallel. (3) We then perform 
a site-resolved measurement of the onsite particle num- 
ber n^ c in each copy, modulo n. We can then deter- 
mine the measured value of the swap operator for 
all possible subsystems being swapped, in parallel 
(where we note that 7Z can also denote the whole sys- 
tem), as these commute. The possible measurement out- 
comes for the swap operations V™, = 1 • • • n} 
can then be computed from the number measurements 
as Yljen e1 ^ nj ' c ^ c ~ 1 \ We illustrate this scheme for 
n = 2 below, and n = 3 in the supplementary material. 

Scheme for n = 2 — As an example, the measurement 
of the Renyi entropy for n = 2 [22 is illustrated in Fig.l. 
In step (2), the FT for two copies is a beamsplitter oper- 
ation, 

ai,i -> + a i>2 )/>/2, a ij2 (a i)2 - a i:1 )/V2. (2) 

This can be achieved by lowering the barrier between 
the two copies (e.g., using a superlattice [35 ) and al- 
lowing the atoms to tunnel from one site to its copy for 
a time T = 7r/(4Ji2) where J12 is the tunnelling rate. 
For this step to work in this form, interactions between 
atoms have to be turned off during this process, e.g. via 
a Feshbach resonance. We will show below that this re- 
quirement can be relaxed in an alternative method for 
hard-core bosons. The number measurement modulo 2 in 
step (3) is then a parity measurement of the site-resolved 
occupation number in each site and copy {n^ c }, exactly 
as is performed in recent quantum gas microscope exper- 
iments mnj. 

The measured value of the swap operator, V^, can 
be computed as (— l)^e^ ni - 2 ; i.e., simply by determin- 
ing whether the total atom number in block 7Z of copy 2 
after the measurement is even or odd. This process is 
illustrated in Fig. lb for one measurement instance, and 
repeating this process allows the expectation value (V^ 1 ) 
to be computed. To obtain this relationship between the 
occupation numbers after the beamsplitter operation and 
(V2 1 ), we note that each possible measurement outcome 
{rii yC } corresponds to the measurement of one of the two 
eigenvalues of the (hermitean) operators V^, which are 
±1. This can be seen as follows: The eigenspaces of 
are the subspaces of the total Hilbert space that are (anti- 
symmetric with respect to exchange of the two copies, 
which we denote by Ti^ (H-^). Since the swap on a sub- 
system 7Z can be constructed by local swaps of the sites 
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FIG. 2: Entanglement buildup and measurement, both for (i) 
softcore bosons in ID after a quench from a Mott Insulator at 
JJ I J — 10 to a superfluid at U/ J = 1; and (ii) hardcore bosons 
with tunnelling J, beginning from an initial state where every 
second lattice site is occupied, (a) Quenches in small systems 
with 8 particles on 8 lattice sites, showing Von Neumann and 
Renyi entropies for a bipartite splitting in the center of the 
system; (b) the same as (a), but for a system of 30 particles 
on 30 sites, and also including 2S2 — S3 (dotted line); (c) The 
number of single shot measurements (#) required to deter- 
mine S n for the Bose-Hubbard quench in (b) with a relative 
accuracy a. For 2S2 — S3, measurements are distributed be- 
tween (S2) and (S3) to minimize the total number, (d) Limi- 
tations of t-DMRG simulations for the Bose-Hubbard quench 
based on different bond dimensions D [10], shown here with 
D = 2 l for 4 < / < 9, with corresponding bound Svn = I.. 
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only a single site. Now, if we denote the annihilation op- 
erator for bosons on site i of copy c by a^ c , c G {1,2}, 
the tunnelling procedure in step (2) gives us Eq. pi). 
This maps the symmetric subspace of the two modes 
and cij y 2 to the subspace of states with an even number 
of atoms in mode 0^2- The anti-symmetric subspace is 
mapped to states with an odd number of bosons in mode 
<2j 5 2- This can be seen by noting that the anti-symmetric 
subspace 7-LJ of the two modes at site j is spanned by the 



states {(at x — at 2 ) 2n+1 |vac)}, while the symmetric one 
is spanned by {(at x -a] ?2 ) 2n |vac), (at x + a] ?2 ) n |vac)} 
(where n = 0,1,2,...). Under the above operation the 
first set is mapped onto {(a^ 2 ) 2n+1 |vac)}, i.e. states 
with an odd number of atoms in a^, while the second 
is mapped onto {(a^ 2 ) 2n |vac), (a^ 1 ) n |vac)} with an even 
number of atoms in aj^- Therefore, the measurement 
outcome of is +1 if nj^ is even and —1 if 77,^2 is 
odd. 

Example of a Mott Insulator- Superfluid Quench — We 
now illustrate the entropy measurement for a quench 
in the Bose-Hubbard model, where quantum dynamics 



generates substantial entanglement in a short time [36]- 
[38] . In Fig. 2 we plot results obtained via t-DMRG 
calculations from dynamics (in each copy of the sys- 
tem) described by the Bose-Hubbard model (h = 1), 

H B h = -JJ2(ij) a \ a 3 + ( u / 2 ) J2i a T a2 n where J is the 
tunnelling rate between neighbouring sites (...), and U 
is the onsite interaction strength. We plot both (i) a 
quench for softcore bosons from a Mott Insulator state 
with U I J = 10 to a superfluid (SF) at [// J = 1, and (ii) a 
quench for hardcore bosons U —> 00, where we begin from 
an initial state with one particle on every second lattice 
site. Such a state could be produced in a superlattice 
potential [6 , or using recently demonstrated techniques 
to directly remove atoms from an initial Mott Insulator 
state [39] , 

Fig. 2a shows Svn and S2 calculated for a bipartite 
splitting in the center of a small system with 8 parti- 
cles on 8 lattice sites. In this size of system, we see that 
the entanglement entropy saturates as the system ther- 
malizes for softcore bosons [40]. In the hardcore case, 
the system is integrable, and we see large oscillations 
in the entanglement entropies, instead of thermalization 
[4T| [42] . Fig. 2b shows a comparison between the von 
Neumann entropy, Svn, and Renyi entropies S n (for 
n = 2,3,4), after a quench in a larger system with 30 
particles on 30 sites. We see rapid growth in the entan- 
glement of two halves of the system, and we note that the 
Renyi entropies provide relatively good bounds for the 
Von Neumann entropy, especially if we use the the bound 
Svn > 2S2 — S3, which is shown as the dashed line in 
the figure. Fig. 2b shows the number of single shot mea- 
surements required to determine the quantities in Fig. 2b 
with a relative error of a. We note that in these larger 
systems, the t-DMRG simulations we use to compute dy- 
namics are limited to simulating short timescales, since 
the von Neumann entropy increases linearly with time 
[IUCL2]. A matrix product state (MPS) with a bond di- 
mension D [10] is capable of representing a maximum von 
Neumann entropy up to log 2 (D). In Fig. 2c, we show the 
evolution of the entropy as a function of time for D = 2 l 
with 4 < I < 9, showing clearly the different bounds for 
the entanglement that can be represented. In the case of 
D = 512 we can faithfully simulate time-scales only up 
to tJ ~ 3. The simulation times scale ~ D 3 , and in an 
experiment, substantially higher entanglement entropies 
could be generated than are accessible in reasonable time 
on a classical computer. This could be demonstrated di- 
rectly by measuring S n (p) in an experiment. 

Simplified scheme for hard-core bosons — The measure- 
ment scheme presented above relies on the ability to turn 
off interactions between the atoms in order to realize 
the beamsplitter operations, which can be challenging 
for some atomic species. However, in certain cases this 
requirement can be relaxed, e.g., in the case of hardcore 
bosons U ^> J, where we have at most one atom per 
site, the measurement can be performed without switch- 
ing the interaction strength. For measurement of S^p), 
the symmetric subspace at site j, is spanned by 



4 



5% 
A 
2.5% 

0% 



tJ = 1.5 




V — tj = i.o 








\^^U=^5 


^ff=o 



5% 
A 
2.5% 

0% 





1 (h\ 


*J= 1.5 




£J = 1.0 - 




= 0.5^ 









-0.02 -0.01 0.01 0.02 0.02 0.04 0.06 0.08 0.1 
e u e 



°a...Af) 



FIG. 3: Errors A = (tr(p 2 a ... M ) - (Uf =a e^)t) /tr(pg 
that are introduced by an imperfect beamsplitter operation, 
determined via t-DMRG simulation of the Bose-Hubbard 
quench in Fig. 2a and resulting measurement process, (a) 
the effect of timing errors, with T — 7r/(4Ji2) + e; (b) error 
introduced by a finite interaction strength u e — U e / J 12 during 
the measurement. 




FIG. 4: Deviations resulting from slightly 

different evolutions in the two copies, A = 

(tr(pl.. M )-(Uf = a^ Aj h)/tT(pL. M )- (a) Results for 
slightly different tunneling amplitudes in each copy after the 
MI-SF quench from Fig. 2a. (b) Results for a case where each 
copy has an additional harmonic trapping potential with 
LJi = 2 x 10J(iV/2) -2 for the first copy, and we otherwise 
perform the same quench as in (a). 



{|vac), (at 1 + at 2 )|vac), at x at 2 |vac)} while the antisym- 
metric one, HJ , is spanned by (at }1 — at j2 )|vac). Of those 
four states the tunnel coupling in step (2) of the measure- 
ment scheme only affects states with one particle in total 
on the copies of the lattice site. The other two states are 
invariant (either because there are no atoms or because 
tunneling is suppressed due to the hardcore constraint). 
Thus we map — » {|vac), at -J vac), at^at^lvac)} and 

V~ — >• at 2 1 vac) in step (2). Measurement of the onsite 
atom number can then be used to directly distinguish be- 
tween the two eigenspaces of and thus determine 
the measurement result. We extend this to the measure- 
ment of Ss(p) in the supplementary material. 

Robustness against imperfections — We now analyze 
the robustness of the measurements to errors in the mea- 
surement steps, especially imperfect implementation of 
the beamsplitter operations in step (2). Two main im- 
perfections can occur here, arising from timing errors or 
residual interactions. In the case of timing errors where 
T = 7r/(4Ji2) + e, we can show analytically [25] that the 
measured value of entanglement will always be smaller 



than the actual entanglement by an amount proportional 
to e 2 , thus leaving a clear lower bound. In Fig. 3a we plot 
results obtained from a full t-DMRG simulation of the 
Bose Hubbard quench from Fig. 2, and the subsequent 
measurement operation. For timing errors on the order of 
1% we find a A on the order of 1% which increases slowly 
with increasing entanglement. In the same calculation, 
we also consider errors introduced by residual interparti- 
cle interactions U e /J\2 present during the measurement 
operation. From Fig. 3b we see that even in the case of 
an interaction of 10% of the beam-splitter tunneling am- 
plitude J 12 we find resulting errors on the order of only 
1% which increase slowly with time. 

Non-identical copies and verification of a quantum 
simulator — Another type of imperfection, is where the 
"copies" undergo different dynamics before the measure- 
ment process, e.g., due to different Hamiltonians. Such 
errors are relevant for characterising the accuracy of the 
quantum simulator itself, as well as constituting an error 
in the measurement scheme. In Fig. 4a we show results 
when we have slightly different tunneling amplitudes in 
each copy after the Bose-Hubbard quench. We find that 
the error introduced by 2% deviations in J correspond 
to a A of 1% on a time-scale of t J ~ 1, which increases 
with time. Fig. 4b shows results from a trapped case, 
where we begin with slightly different trapping poten- 
tials in the individual copies, and then perform the same 
interaction quench from U/J = 10 to U/J = 1. On a 
time-scale t J ~ 1 errors of 1% lead to deviations on the 
order of 0.1%. In this case, the quantity that is actu- 
ally measured is tr{pip2}, where p\ and p2 are the den- 
sity operators of the different copies. This can be used 
as a tool to verify the quantum simulator in the case 
that the evolutions should be identical, or to investigate 
dynamics with the evolution Hamiltonians being set to 
slightly different values. The latter could be interesting 
in regimes of the Bose-Hubbard model expected to dis- 
play signatures of quantum chaos [14, 43 , 44 , where the 
inner product might decay exponentially in analogy with 
single-particle systems exhibiting quantum chaos [45HI7] . 

Summary — We have analyzed measurement of Renyi 
entropies of ID bosons in an optical lattice during a 
quench, using techniques currently available in quantum 
gas microscope experiments. Such a tool can be used to 
verify the coherence and accuracy of a ID and 2D quan- 
tum simulator, and determine the entanglement growth 
in a quantum quench, helping the experiments to diag- 
nose regimes where dynamics can be realized that are not 
accessible to present classical computations. 
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Supplementary material 

I. RELATIONSHIP BETWEEN THE RENYI 
ENTROPY AND THE EXPECTATION VALUE 
OF THE SWAP OPERATOR 

Let us assume that we have two states p\ and p2 
in separated subsystems 1 and 2 with identical Hilbert 
spaces (these represent the two copies in our measure- 
ment scheme). We then define a swap gate V2 on these 
two states, so that for any two states of the systems |^i)i 
and ^2)2, 

Vs#i)i <8> ^2)2 = ^2)1 <8> ^1)2. 

We then decompose the density operators of the subsys- 
tems as pi = Y.ijPifl^Ul P2 = T,kiPki\ k )( l \ in the 
corresponding basis states for each subsystem. We then 
see that 

tr{V 2P(1) ®p (2) } = tr\v 2 Y l P$ ) P ( ?MV\ ® \ k )( l \ I 

[ ijkl ) 

[ijkl J 

ijkl 

= ^PikPki = tr{pip 2 }- 

ik 

If p = pi = p2 (i.e., we have two identical copies), this 
results in tr{V2/0(8> p} = tr{p 2 }. This result can be read- 
ily generalized to the case of permuation of n copies, to 
obtain tr{p n } = tr{V n p® n }, where the shift operator on 
n copies, V^i) . . . \ijj n ) = \^ n )\^i) • • • |^n-i)- 



II. REQUIRED NUMBER OF 
MEASUREMENTS TO DETERMINE S n (p) WITH 
PRESCRIBED ACCURACY 

The number of measurements needed to determine the 
expectation value of with a given accuracy can be 
determined by calculating its variance, which can be cal- 
culated in terms of (V™) by noting that V^V^ = 1. For 
n > 2 the shift operator V n is no longer hermitian. There- 
fore in principle it is not guaranteed that its expectation 
value is real. However, for our assumption, that we have 
a state of n copies (p n ), it is. Assuming that this require- 
ment is fulfilled, we need to determine only the real part 
of its expectation value. (In other words, measure the ex- 
pectation value of (V n -\-V 7 l)/2). This reduces the number 
of measurements, by roughly a factor of 1/2. (However, 
if one determines also the imaginary part of (V n ), this as- 
sumption can be tested.) The number of measurements 
needed to determine (V 7 ^) with a relative error a is then 



given by = (1/2 - (V^) 2 /2)/(a(V^)) 2 if n is odd, 

and = (1/2 + <0 2 /2 " < if « » 

even. This shows that the effort to determine Renyi en- 
tropies of higher orders grows exponentially with n (Note 
that one has tr{p m } ^ (tr{p n }) m / n if n > ra). However, 
measurement of S n (p) for several low values of n can al- 
ready be useful in finding bounds for the von Neumann 
entropy Svn- 



III. MEASUREMENT OF THE n = 3 RENYI 
ENTROPY 

Here we give the specific example for implementation 
of the quantum Fourier transform (FT) for measure- 
ment of the n = 3 Renyi entropy. The unitary map 
associated with the FT is generated in three steps, 
that involve tunnelling between the copies, and shifting 
the relative potential depths between copy p and v\ 
T^ v (a, 4>) : —> a M cos(a) — a v sin(a); a v — >• (a v cos(a) + 
a M sin(a))e _z ^;. One has (up to local phases) U% T = 

T 2 , 1 (^/4,V2)T 3 ,i(arctan(l/v / 2),2V3)T 3 ,2(-V4^/3). 
Note that the second step, the "beam splitter" between 
copy 1 and 3 (which are not neighboring) can be broken 
down to three beam splitting operations that involve 
only neighboring copies. 

For arbitrary n the FT can be achieved with n{n — l)/2 
beam splitter operations [34 j. If only beam splitters 
between neighboring sites are allowed one can get it 
with [n — l)n(2n — l)/6 of these operations. 



IV. SIMPLIFIED MEASUREMENT OF S 3 (p) 
FOR HARDCORE BOSONS 

Here we show that we can also measure tr{p 3 } for 
hardcore bosons without switching off the interparti- 
cle interactions. To do so, we apply the same se- 
quence of lattice modulations that would produce the 
U% T in the absence of interactions. The states |vac) 
and a]^a]^ a ],3 1 vac ) (correspondingto the eigenspace of 
with eigenvalue one), are unaffected by this. 
The eigenstates states with one particle in total on 
the three sites, ^= Y^i=i l e z ^^ l ~ 1 ^ k ~ 1 " > |vac) for k = 
1,2,3 (which are eigenstates with eigenvalue e^^ -1 )) 
are mapped onto the three states at fe |vac). This is 
the same as in the general scheme, as there is just 
one particle, which consequently does not feel the 
hardcore constraint. From the particle- hole symme- 
try we deduce that the remaining three basis states 

(with eigenvalue e*^^ -1 )), which have one hole, are 
mapped into the states that have the hole in copy fc, 
a j,k a j i a ] 2 a \ 3l vac )- Clearly all these states can be dis- 
tinguished in a measurement of the local occupation num- 
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ber. Note that, in contrast to the above general scheme The measured expectation value is then 
where one requires the ability to measure the onsite parti- 
cle number modulo three, this needs only a measurement 
modulo two. 



V. ANALYSIS OF TIMING ERRORS IN THE 
BEAMSPLITTER OPERATION 



tr {V?{e)p ® p) « tr {p 2 n } + e 2 tr {<D n Vfp p) (5) 



Here we analyse errors in the beamsplitter operation in 
step (2) of the measurement process. We assume that the 
"rotation angle" a is not perfectly chosen: a = 7r/4 + e, 
e«l, resulting in a operation 



1 



a j,l + a j,2 a j,l ~ a j,2 



V2 
V2 



72 ■ 



(3) 
(4) 



Note that, since V^p p is a positive operator and 
O n = 2 ^j G 7^( a ],i a i,2 - a ],2 a i,i)) is a negative one, 

these errors always lead to a measured value of tr {/^} 
that is smaller than the actual one. 



